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I.  INTRODUCTION 


A.  Background 

Since  the  mid  1950's  and  especially  the  last  five  years,  considerable 
effort  has  been  expended  by  various  DOD  agencies,  the  Ballistic  Research 
Laboratory  being  one,  to  measure  the  motion  of  gun  tubes  during  the  interior 
ballistic  and  launch  cycle.  The  primary  purpose  of  these  measurements  is  to 
obtain  correlations  between  gun  system  dynamic  behavior  and  overall  system 
accuracy.  The  secondary  purpose  is  to  validate  numerical  gun  system 
simulations  to  predict  system  behavior.  Although  various  measurement 
techniques  have  been  tried  with  varying  degrees  of  success,  the  use  of 
accelerometers  still  holds  the  most  promise,  provided  their  measurement  results 
can  be  unequivocally  interpreted.  The  primary  advantage  of  accelerometers  is 

that  they  can  be  applied  under  field  conditions  with  limited  restriction  or 
modification  of  the  gun  system. 

B.  The  Basic  Problem 


The  basic  problem  is  to  measure  the  components  of  dynamic  motion  of  the 
gun  tube  that  significantly  affect  the  launch  conditions  of  the  projectile  and 
its  subsequent  terminal  accuracy.  The  measurement  problem  becomes  one  of 
being  able  to  discern  the  significant  vector  components  of  motion  at  a  point 
on  the  gun  tube  from  a  complex  combination  of  vector  components  generated  by 
various  dynamic  phenomena,  including  the  ones  of  interest.  Unlike  well- 
controlled  laboratory  experiments  designed  to  separate  various  physical 
phenomena,  gun  system  acceleration  measurements  are  subject  to  interference 
from  local  accelerations  generated  by  stress  waves,  dilational  vibrations, 
traveling  loads,  impulsive  forces  and  unbalanced  projectiles. 


c«  Approach  to  the  Problem 

To  gain  a  meaningful  measurement  from  accelerometers,  one  must  take 
advantage  of  the  characteristics  of  the  physical  phenomena  encountered  during 
the  measurement.  Remembering  that  the  global  structural  motion  is  paramount, 
emphasis  is  placed  on  the  low-frequency  response  of  the  system.  By  separating 
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analytically  and  by  appropriate  filtering,  the  high-frequency  interference  can 
be  eliminated  from  the  measurement.  Separation  of  the  desired  vector 
components  of  motion  is  achieved  through  the  judicious  design  of  sensor  arrays 
based  on  vector  mechanics. 

1.  Fourier  Analysis  and  Cross  Correlation  Technique.  These  techniques 
are  applied  during  the  data  analysis  and  rely  on  the  spectral  nature  of  the 
physical  phenomena  being  measured  to  facilitate  separation  of  effects •  The 
viability  of  these  techniques  is  strengthened  by  proper  sensor  array  design 
and  filtering  techniques . 

2.  Mechanical  and  Electronic  Filtering  Techniques.  By  mechanical 

filtering  between  the  sensor  and  the  phenomena  and  electronic  filtering 

between  the  sensor  and  the  acquisition  system,  undesirable  high  frequency 

physical  effects  can  be  separated  from  the  measurement.  This  allows  greater 

sensitivity  to  be  achieved  in  the  measurement  with  increased  sensor 

reliability  and  decreased  burden  on  the  analysis.  However,  these  techniques 

must  be  carefully  designed  for  the  specific  phenomena  to  be  measured  and  for 

2  3 

the  specific  system  on  which  the  measurement  is  made.  y 

3.  Spacial  Arrays  and  Vector  Discrimination  Techniques.  To  implement  a 
viable  acceleration  measurement  in  a  six-degree-of-f reedom  environment, 
sensors  must  be  mounted  in  a  spacial  array  in  order  to  discriminate  and 
separate  the  desired  vector  components  from  the  measured  data.  The 


i 


^James  N.  Walbert,  "Application  of  Digital  Filters  and 
Fourier  Transform  to  the  Analysis  of  Ballistic  Data, " 
Technical  Report  ARBRL-TR-02347 ,  1981  (ADA  102890). 


the 

BRL 


^ James  N.  Walbert,  "Computer  Algorithms  for  the  Design  and. 
Implementation  of  Linear  Phase  Finite  Impulse  Response  Digital 
Filters ,  "  BRL  Technical  Report  ARBRL-TR-0234  6 ,  1981  (ADA  103112). 


^  James  0.  Pilcher  II,  "Application  of  Mechanical  Filters  to 
Ballistic  Measurements ,  "  BRL  Technical  Report  (to  be  published) 
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application  of  these  techniques  governs  the  required  analyses  techniques  and 
enhances  the  efficiency  of  filtering  techniques. 

The  remainder  of  this  report  discusses  the  analytical  rationale  behind  the 
design  of  spacial  arrays  to  facilitate  separation  of  the  vector  components  of 
gun  tube  motion.  The  other  techniques  mentioned  above  will  be  discussed  in 
subsequent  reports. 

II.  SIX-DEGREE-OF-FREEDOM  MOTION  OF  A  POINT  ON  A  GUN  TUBE 

In  order  to  design  an  accelerometer  array  for  measuring  the  motion  of  a 
gun  tube,  a  description  of  the  motion  of  the  point  of  measurement  is 
required.  Once  such  a  description  is  achieved,  arrays  can  be  designed  to 
enhance  or  minimize  certain  vector  components  of  acceleration  as  desired. 
Naturally,  any  such  description  must  rely  on  assumptions  about  the  behavior  of 
the  gun  tube. 

A.  Assumptions 


The  following  are  the  basic  assumptions  about  the  behavior  of  a  gun  tube 
and  the  resulting  motion  at  a  point  on  the  surface  of  the  gun  tube. 

1*  Cross- Sectional  Property .  A  plane  cross  section  of  the  gun  tube 
remains  planar  throughout  the  duration  of  the  dynamic  response  of  the  gun 
tube.  This  assumes  that  local  bending  deformations  are  very  small  and 
elastic.  Phenomena  that  violate  this  assumption  are  passage  of  a  cocked 
projectile,  load  or  geometric  discontinuities  in  the  immediate  vicinity  of  the 
point  of  measurement  and  unbalanced  projectiles.  Figure  1  illustrates  this 
assumption • 

2*  In-Plane  Behavior  of  The  Cross  Section.  The  cross  section  of  the  gun 
tube  can  undergo  both  symmetrical  and  asymmetrical  deformation  about  the 
central  axis  in  the  plane  of  the  cross  section.  This  assumption  allows 
assessment  of  the  acceleration  components  generated  by  dilational  vibration, 
ovality  due  to  local  bending.  Bourdon  effects  and  projectile  passage.  Figure 
2  illustrates  this  assumption. 
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Figure  1.  Planar  Properties  of  a  Gun  Tube 


SYMMETRIC  BEHAVIOR 


ASYMMETRIC  BEHAVIOR 


Figure  2.  In-Plane  Behavior  of  a  Gun  Tube 
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Local  Coordinate  System.  The  local  coordinate  system  is  assumed  to 
be  Cartesian, fixed  with  respect  to  the  cross-section  plane,  with  two  of  the 
axes  lying  in  the  plane  and  the  origin  at  the  center  line  of  the  tube.  This 
coordinate  system  translates  and  rotates  with  respect  to  a  fixed  Earth 
coordinate  system.  Consequently,  any  vector  that  is  directionally  fixed  in 
the  plane  of  the  cross  section  has  a  time  derivative  of  the  unit  vectors  due 
to  their  changes  in  direction.  Figure  3  illustrates  this  assumption. 


X 


z 


Figure  3.  Local  Coordinate  System 

^ ’  Continuity .  The  motion  is  assumed  to  have  a  continuous  second  time 
derivative.  That  is,  the  third  time  derivative  is  defined  and  bounded  except 
at  discrete  points  where  it  possesses  left-  and  right-hand  limits.  This 
assumption  is  based  on  the  gun  structure  conforming  to  the  laws  of  classical 
mechanics  of  rigid  and  deformable  bodies.  Deviations  from  this  assumption 
would  be  indicative  of  failure  of  a  structural  component  or  the  sensor. 


Inertial  Reference.  Since  accelerometers  are  inertial  devices,  the 
measurement  constitutes  an  observation  from  a  point  displaced  away  from  the 
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center  of  rotation  of  a  body  moving  through  space.  The  sensed  accelerations, 
then,  are  the  components  of  the  translational  motion  dependent  on  the 
instantaneous  angle  of  rotation  of  the  local  coordinate  system  with  respect  to 
the  Earth  coordinate  system,  in  conjunction  with  the  cross  products  of  the 
local  displacement  vector  with  the  angular  rates  of  change  of  the  local 
coordinate  system. ^ 

B.  Position  Vector 


The  position  of  the  point  of  measurement  on  the  outer  surface  of  the  gun 
tube  is  described  in  terms  of  the  unit  vectors  Ux,  Uy  and  Uz  of  the  local 
coordinate  system  by  the  position  vector, 

P  =  R  +  7  ,  (1) 

where  R  =  the  radius  vector  emanating  from  the  origin  of  the  Earth 

coordinate  system  and  terminating  at  the  origin  of  the  local 
coordinate  system. 

r  =  the  radius  vector  emanating  from  the  origin  of  the  local 

coordinate  system  and  terminating  at  the  point  of  measurement. 


Figure  4.  Position  Vector 


S.W.  MoCuskey,  IntmdMeti^n  to  Advanced  Dynamais,  Addison-Wesley  Publishing 
Company,  1958. 
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In  expanded  form, 


p  =  <V“rx>  ux  +  ^VV  uy  +  (Rz+rz)  uz-  (2) 

C.  Velocity  Vector 

The  velocity  at  the  point  of  measurement  on  the  tube  is  obtained  by 
forming  the  total  time  derivative  of  the  position  vector.  Bear  in  mind  that 
the  derivative  of  the  unit  vector  U  exists  because  of  the  rotation  of  the 
local  coordinate  system  and  is 


dU 

dt 


a)  x  U 


where  a>  =  the  angular  rate  of  change  of  U  which  is  of  constant  magnitude 


Therefore , 


V 


dP 

dt 


=  R  +  r+  u)x  (R  +  r) 


(3) 


where 


a)  =  the  time  rate  of  change  of  the  local  coordinate  system; 


a)  x  r  =  is  the  velocity  normal  to  the  position  vector,  r. 


In  expanded  form, 
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V 


=  [R  +r  +  (R  +r  )io  -  (R  +r  )a>  ]U 
xx  z  z  y  y  y  z  x 


+  [R  +r  +  (R  +r  )u>  -  (R  +r  )ui  ]U 

yy  xxz  zzx  y 


+  [R  +r  +(R  +r  )w  -  (R  +r  )u>  ]U 
z  z  y  y  x  x  x  y  z 


(4) 


The  term  r  in  Eq .  (3)  results  from  the  flexibility  of  the  tube  cross  section 
and  vanishes  for  a  rigid  cross  section.  Figure  5  illustrates  the  velocity 
vector,  V. 


Figure  5.  Velocity  Vector 


D.  Acceleration  Vector 


The  acceleration  vector  is  the  total  time  derivative  of  the  velocity 
vector,  V. 
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d_V  _  dR  ,  d_r  ,  d((^x(R+r)  ) 
dt  dt  dt  dt 


•  •  ••  " 

A  =  R  +  r  +  2 ( o)x ( R+r )  4*  a)(w  (R+r ) )-(R+r ) ( o>*o>)  +  ojx(Rfr)  . 

(5) 


where 


2 ( cox ( R+r ) )  +  (A)(u>*R+r)  -  (R+r)(u>*u))  is  the  coriolis  accelerations; 


u>x(R+r)  is  the  tangential  accelerations. 


In  expanded  form. 


A 


) 


99 

+  [R+r  +o)  a)  (R+r  )  +o)  a)  (R+r  )-(R+r  )(u>  +u)  ) 

y  y  y  z  z  z  y  x  x  xy  v  y  y/v  z  x  ' 


+  2o)  (R+r  )-2o)  (R+r  )+a)  (R+r  )-o)  (R  +r)]U 
z  X  X  x  z  z7  zv  X  X  X  z  2J  V 


+ 


2_^  2 
(A)  +0) 

X  y 


) 


(6) 


The  expanded  acceleration  equation,  Eq.  (6),  demonstrates  the  problem  of 
separating  specific  vector  components For  studying  and  determining  accuracy 
effects,  the  vector  components,  u>,  R,  R,  in  terms  of  the  Earth  coordinate 
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system,  are  generally  the  primary  components  of  interest .  These  vectors  are 
derived  by  attempting  to  separate  and  transform  R  and  03  and  evaluating  their 
subsequent  integrals.  These  components  are  often  used  to  verify  numerical 
simulation  results.  Rather  than  follow  these  techniques,  a  different  approach 
is  in  order. 

The  angular  velocity  square  terms  and  the  inner  product  of  angular 
velocity  constitute  the  products  of  frequency  terras  which  generate  beats  and 
second  harmonics.  This  means  that  measurement  and  acquisition  system 
bandwidths  and  data  rates  must  be  commensurate  for  frequency  response  twice  as 
high  as  those  predicted  by  normal  mode  analyses  • 

E.  Some  Simplifications 


Before  progressing  further  an  examination  of  the  terms  in  Eq .  (6)  is  in 


order.  By  rearranging  Eq .  (6)  and  substituting 


M 


x 


[R  +03  0)  R  +o>  03  R  -R 
x  xyy  xzz  x 


2 


)+2o3  R  -2u)  R 
y  z  z  y 


+0)  R  -o)  R  ]  ; 
y  z  z  y 


M 


[R  +0)  03  R  +03  03  R-R 

y  yzz  yxx  y 


(w 


2 


)+2o3  R  -2o3  R  +o)  R  -03  R 
zx  xzzxxz 


> 


M  = 
z 


[R  +03  0)  R  +0)  0)  R  -R 
z  zxx  zyy  z 


(w 


)  +  203  R  -2o3  R  +03  R  -0)  R  ] 
xy  y  x  x  y  yx 


A 


**  2  2 

[M  +r  +o)  03  r  +o)  03  r  -r  (o)  +U3  )+2o3  r  -2o)  r  +03  r  -o>  r  1 U 

xxxyyxzzxy  z  yz  zyyzzy  x 


+ 


[M  +r  +o)  oj  r  +o)  o)  r  -r  (03 
y  y  yzz  yxx  y  z 


)+2o)  r  -2 u»  r  +u>  r  -u  r  ]U 
zx  xz  zx  xz  y 


+  [M+r+oiojr+ojojr 
z  z  zxx  zyy 


22v  •  ••  •  — 

( 03  +03  )+2o3  r  -2o3  r  +03  r  -03  r  ]U 

zx  y  xy  y  x  x  y  y  x  z 


(7) 
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The  terms  My,  Mz  are  the  components  of  relative  motion  between  the  origins 
of  the  Earth  and  local  coordinate  systems. 


Eq.  (7)  can  be  immediately  simplified  by  applying  the  previous 
assumptions.  Since  the  measuring  point  lies  in  the  same  plane  as  the  origin 
and  two  axes  of  the  local  coordinate  system,  its  description  becomes  two- 
dimensional.  As  a  matter  of  convention,  the  coordinate  axes  often  used  are 
the  x  axis  in  the  vertical  direction,  the  y  axis  in  the  horizontal  direction 
and  the  z  axis  perpendicular  to  the  x-y  plane  and  tangent  to  the  axis  of  the 
gun  tube  at  its  intersection  with  the  x-y  plane.  This  means  that  the  tube 
cross  section  at  the  point  of  measurement  lies  in  the  x-y  plane  and  that  the 
rz  component  of  the  position  vector  has  zero  magnitude.  This  immediately 
simplifies  Eq.  (7)  to  that  shown  in  Eq.  (8). 


(8) 


When  the  sensor  is  placed  on  either  the  x  or  y  axis,  Eq.  (7)  reduces  to  Eqs . 


(9)  and  (10). 


When  the  sensor  is  on  the  y  axis,  r  ,  r  and  r  are  equal  to  0; 

X  X 


thus 


(9) 


When  the  sensor  is  on  the  x  axis,  ry,  ry  and  ry  equal  0;  thus, 
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••  ••  90 

=  [M  +r  -r  (w  ~Hi)  )  ]  U 
x  x  x  y  z  z 


+  [M  +o)  a)  r  +2o)  r  +0)  r  ]  U 
y  yxx  zx  zx  y 


•  •  — 

4-  [M  +0)  a)  r  -a)  r  -a)  r  )U  . 
z  zxx  yx  yx  z 


(10) 


Further  simplification  can  be  achieved  by  taking  advantage  of  the  character  of 
the  position  vector  r.  Two  types  of  motion  occur  in  the  cross  section  of  the 
tube,  symmetrical  and  asymmetrical.  Both  of  these  types  of  motion  are  radial 
deflections  induced  by  local  strain  effects  that  tend  to  be  much  higher  in 
frequency  than  the  structural  response.  Though  these  displacements  are  small 
enough  to  ignore,  their  corresponding  accelerations  are  quite  intense.  These 
phenomena,  with  the  exception  of  the  impulse  generated  by  projectile  passage 
through  the  plane,  generate  high-frequency  and  high-intensity  oscillating 
accelerations.  These  oscillations  are  in  excess  of  2  kg's  in  amplitude  and 
above  7  kHz  in  frequency.  Consequently,  these  oscillations  can  be  eliminated 
from  the  data  by  mechanical  filtering  between  the  sensor  and  the  gun  tube. 

The  vector  r  can  be  recast  as 

r  =  d  +  S  +  £ ,  (11) 

where 


d  =  gun  tube  radius  in  plane  of  measurement  and  is  of  constant  length; 


6  =  oscillatory  displacement  about  the  length  of  d; 


£  =  impulsive  displacement  about  d. 
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The  impulse  displacement  £  can  often  be  separated  by  Fourier  techniques  during 
the  data  analyses.  Thus,  by  eliminating  the  effects  6  and  C  from  vector  r, 
Eqs.  (9)  and  (10)  reduce  to  Eqs.  (12), 


For  r  =  r  U  , 
x  x 


_  2  2”  «  _ 

A  =  [M  -r  (u>  +«»  )]U  +[M  -Ho  a)  r  -Ho  r  ]U  +[M  -Ho  id  r  -u>  r  ]U  ; 

XXY  z  x  yyxxzx  yLzzxxyx  z* 


for  r  =  r  U  , 

y  y 


A 


[M  +oj  a)  r  -a)  r  ]U  +[M  -r  (oj  ^)] 
x  xyy  z  y  x  y  y  y  zyj 


U  +[M  -fa)  oj  r  -hL  r 
y  z  z  y  y  x  y 


(12) 


Further  simplification  of  the  problem  involves  sensor  selection  and 
directional  properties • 

III,  GENERAL  EQUATION  FOR  ACCELEROMETER  RESPONSE 

Accelerometers  are  not  unidirectional5  devices;  their  response  is  three 
dimensional  in  character.  Because  most  accelerometer  measurements  are 
performed  under  conditions  of  relatively  unidirectional  excitation  such  as 
vibration  and  shock  machines,  the  assumption  of  unidirectionality  is  often 
made  without  appreciably  degrading  the  measurement.  However,  under  the 
conditions  arising  from  gun  and  projectile  dynamics,  such  an  assumption  is 
very  misleading  and  can  cause  considerable  degradation  of  the  measurement.  In 
fact,  most  measurements  cannot  be  made  without  taking  advantage  of  the  three- 


Charles  C .  Crede  and  CyviZ  M .  Harris ^  Shock  and  Vibration  Handbook ,  Volume 
ly  McGraw-Hill,  1961. 
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dimensional  character  of  the  accelerometer.  In  order  to  establish  a  viable 
relationship  for  accelerometer  sensitivity,  the  vector  characteristics  of  the 
sensing  element,  its  misalignments  with  external  reference  surfaces  and  local 
coordinate  systems  must  be  considered. 

A.  Accelerometer  Response  Characteristics 


Assuming  that  the  accelerometer  is  subjected  to  an  environment  within  its 
capabilities  in  both  the  magnitude  and  frequency  domain,  its  response 
characteristics  can  be  expressed  as  a  vector  S(a,T,7r)  where  a,  t,  tt  ,  are  the 
unit  vectors  of  the  sensor  coordinate  system.  This  coordinate  system  is  not 
necessarily  orthogonal,  but  the  deviations  from  orthogonality  are  generally 
very  small  angles.  In  addition,  there  generally  exist  small  rotations  of  the 
principal  components  of  the  sensitivity  vector  with  respect  to  the  external 
reference  system  of  the  device.  The  angles  of  rotation  and  nonorthogonality 
are  sufficiently  small  to  be  lumped  together  in  the  direction  angles 
a,  3,  y.  By  performing  an  orthogonal  transformation , ^  the  components  of 
sensitivity  can  be  expressed  in  terms  of  an  orthogonal  coordinate  system  with 
unit  vectors  x,  y  and  z  which  are  referenced  to  the  external  geometry  of  the 
device.  Figure  6  illustrates  the  relationships  of  these  coordinate  systems. 
Figure  7  shows  a  comparison  of  an  ideal  accelerometer  response  surface  with  a 
realistic  response  surface . 

The  external  sensor  geometry  is  generally  defined  as  having  a  principal 
axis  that  will  be  defined  as  being  colinear  with  z  and  a  mounting  surface 
that  is  a  plane  either  normal  or  parallel  to  the  principal  axis.  In  our 
discussions,  we  will  address  the  normal  plane  case.  Another  condition  placed 
on  the  analysis  is  that  x  lies  in  the  plane  defined  by  tt  and  a.  This  condition 
can  be  readily  established  from  calibration  data.  The  sensitivity  vector  is 
defined  as 


6 


Richmond  B .  McQuistan 3  Scalar  and  Vector  Fields > 


John  Wiley  &  Sons3 


1965. 
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Figure  6.  Geometrical  and  Sensing  Coordinates 


IDERL  RCCELEROMETER 


F^-Sure  7.  Sensitivity  Surfaces  for  Ideal  and  Real  Accelerometers 
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(13) 


S  =  So+TT+Pir=(a  S+a  T+a  P)x+(a  S+a  T+a  P)y 
xa  xt  xtt  y  a  yx  yir 

+(a  S+a  T+a  P)z, 
z  a  ZT  Z7T 


where 


a..  =  is  the  direction  cosine  between  unit  vector  i  and  unit  vector  i. 
ij  J 

Eq.  (13)  is  only  applicable  when  ot,  g,  y  are  sufficiently  small  such  that 
sin(a)=a,  sin(g)=g  and  sin(y)=y  and  cos(a)=l>  cos(g)=l  and  cos(y)=l.  General 
measurement  accuracy  dictates  that  these  angles  should  be  less  than  3.4 
degrees  or  0.06  radian.  Table  1  gives  the  transformation  coefficients. 


TABLE  1.  TRANSFORMATION  COEFFICIENTS  (Aij.) 


_ 

a 

T 

TT 

X 

1 

-g 

a 

y 

g 

1 

-y 

z 

-a 

y 

i 

Eq.  (13)  then  becomes 


S  =  (S-gT+aP)x+(gS+T-yP)y+(-ctS+yT+P)z.  (14) 


Eq •  (14)  can  be  stated  with  simpler  coefficients  relatable  to  the  cross-axis 
sensitivities  as  measured  during  calibration. 


S  =  Pk^x  +  +  Pz* 


(15) 
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where 


P  =  primary  sensitivity  in  units/g; 

=  maximum  cross-axis  sensitivity  in  %/ 100 ; 
k2  =  value  of  cross-axis  sensitivity  orthogonal  to  kj  in  %/100. 

For  some  accelerometers  which  have  nearly  perfect  alignment,  k ^  =  k2 


S  =  Pk^x  +  Pk^y  4-  Pz. 


(16) 


For  an  ideal  accelerometer  k^  and  k2  are  zero;  thus, 


S  =  Pz. 


(17) 


The  conditions  in  Eq.  (16)  are  often  found  by  screening  accelerometers 
with  less  than  1%  cross-axis  sensitivity.  The  conditions  in  Eq.  (17)  have 
been  approached  by  some  piezoresistive  accelerometers,  but  not  without 
degrading  some  other  properties.  Eqs .  (15),  (16)  and  (17)  show  that 
calibrations  describing  the  behavior  in  three  orthogonal  planes  are  required 
for  six— degree— of —freedom  measurements.  The  remainder  of  our  discussion  will 
use  Eq.  (15)  as  the  basis  of  analysis. 

B.  Accelerometer  Mount  Misalignment 

The  misalignment  of  accelerometer  axes  due  to  mounting  conditions  can 
introduce  significant  error  into  the  measurement .  For  accelerometers  that  are 
bonded  or  fastened  to  their  mounting  surface  by  external  means,  mounting  error 
can  be  eliminated  by  the  use  of  precision  fixtures  and  techniques  during  the 
mounting  process.  However,  stud  mounted  accelerometers  limit  the  ability  to 
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align  the  transverse  axis  in  the  desired  direction.  This  type  of  misalignment 
often  causes  large  angles  of  rotation  and  requires  a  different  correction 
technique  than  the  first  type  mentioned. 

1.  Bonded  Accelerometer .  Although  precision  fixtures  and  techniques  are 
used  in  the  mounting  process,  misalignments  do  occur  due  to  human  error. 

Errors  resulting  from  such  misalignments  can  be  corrected  by  the  same  type  of 
transformation  discussed  in  the  previous  section.  Figure  8  shows  the 
accelerometer  coordinates  misaligned  with  the  local  coordinate  system. 


S(x,y,z)  =  S(U1,U2,U3) 


08) 


Figure  8.  Mounting  Misalignment  for  Bonded  Accelerometer 
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Table  2  shows  the  transformation  coefficients  for  the  angles  e,  <5,  E,  shown 
in  Figure  8 . 

TABLE  2.  TRANSFORMATION  COEFFICIENTS  (e,  <5,  E,  ) 


x  y  z 


S  =  PCkj-e^-SJUj+PCekj+k^fi^+p  (^-6^+1  )U3 , 


(19) 


or, 


S  =  al^+bi^+cl^, 


(20) 


where 


a  =  P(k1-  ek2-  £); 
b  =  P(  ek^+  k2+  6  ) ; 
c  =  P(5k1-  6k2+l). 
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2.  Stud  Mounted  Accelerometers.  Stud  mounted  accelerometers  pose  the 
problem  that  the  final  alignment  of  the  accelerometer  may  include  a  large 
angle  of  misalignment  which  is  usually  in  the  cross-axis  plane.  Because  of 
the  large  angle,  a  two-step  transformation  is  required.  First,  an  arbitrary 
local  coordinate  system  is  assumed  with  one  of  its  unit  vectors  colinear  with 
the  maximum  cross-axis  sensitivity  vector  as  shown  in  Figure  9.  An  orthogonal 
transformation  is  performed  resolving  S(x,y,z)  into  SCU'.U'.U').  Table  3 
gives  the  coefficients  of  this  transformation. 


TABLE  3.  TRANSFORMATION  COEFFICIENTS  (e=0) 


X 

y 

z 

V 

i 

0 

-5 

V 

0 

l 

6 

V 

5 

-6 

1 

Figure  9.  Mounting  Misalignment  for  Stud  Mounted  Accelerometer  to  Adjusted 
Local  Coordinates 
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The  transformed  vector  is 


s  =  P(k1-5)u1+p(k2+6)u2  +  PCSkj-fik  +i)u3. 


(21) 


The  second  step  is  to  perform  the  transformation  of  the  rotation  of  the 
large  angle  e  about  which  is  colinear  with  as  shown  in  Figure  10. 

^1 >  ^2 ,  U3  are  the  local  coordinates.  Since  angle  e  is  large,  cos(e)  ±  1  and 
sin(e)  *  e.  The  coefficients  for  this  transformation  are  given  in  Table  4. 


TABLE  4.  TRANSFORMATION  COEFFICIENTS  FOR  LARGE  e 


V 

V 

V 

*1 

cos (e ) 

-sin(e) 

0 

Uo 

4. 

sin(e) 

cos(e) 

0 

C3 

0 

0 

1 

i  L 


Figure  10.  Mounting  Misalignment  for  Stud  Mounted  Accelerometer  to 
Local  Coordinate  System 
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S  =  P[(k1-5)cos(e)-(k2+6)sin(e)]U1 
+  P[(k1-5)sin(e)+(k2+6)cos(e)]U2 

+  P(5k1-6k2  +  1)U3,  (22) 


or, 


S  =  aU1  +  bU2  +  cU3, 


(23) 


where  a  =  P [ (k^-£ )cos ( e )  -  (k2+6 )sin( e ) ]  ; 

b  =  P[ (k  -5)sin(e)  +  (k2+6)cos(e)] ; 
c  =  P(56k  -6k2+l). 


IV.  ACCELEROMETER  ARRAYS 

To  facilitate  separation  of  vector  components,  accelerometers  must  be 
arranged  in  distinct  arrays  about  the  plane  of  measurement  for  a  specific 
situation.  This  allows  the  generation  of  sets  of  vector  equations  required  to 
separate  components  . 

Particular  interest  is  paid  to  the  separation  of  the  relative 
acceleration  term  M^,  and  Mz .  These  terras  are  the  structural  response 

terms  that  are  predicted  by  gun  simulation  programs.  The  terms  involving 
products  of  r,  r  and  r  are  local  effects  generated  by  local  strains  which 
create  accelerations  that  mask  the  raeasurand  and  also  cause  nonlinear  response 
of  the  sensor.  The  strains  caused  by  wave  propagation  are  removed  by 
mechanical  filtering  techiques .  Other  high-frequency  phenomena  can  be  removed 
by  electronic  methods;  however ,  low  frequency  variations  in  r,  particularly  due 
to  projectile  passage,  must  be  separated  by  comparative  measurement  of 
different  positions  in  an  array. 
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A.  Single-Pair  Array 


The  basic  array  is  a  pair  of  accelerometers  displaced  in  space  and 
oriented  to  specific  directions*  Although  the  convenience  of  symmetry  is  used 
in  this  analysis,  it  is  not  necessary  and  can  be  circumvented  by  analysis. 

The  illustration  in  Figure  11  shows  a  pair  of  accelerometers  mounted 
diametrically  opposite  one  another  on  the  surface  of  the  gun  tube.  The 
origins  of  their  coordinate  systems  lie  in  the  plane  of  the  cross  section  with 
their  unit  vectors  parallel  to  the  coordinate  system  of  the  cross  section  and 
antiparallel  to  each  other .  The  displacement  vectors  to  the  center  of  the 
local  coordinates  from  the  center  of  the  tube  cross  section  are  r^  and  r^. 


When 


-  ri^x>  t^le  acce^-erati°n  vector  is 


A 


1 


-a)  r 

y 


-a)  r 

i  y 


z 


(24) 


When 


r3  ~~  r3^x  =  ”  rA>  t*le  acceleration  vector  is 


A  = 


2.  2* , - 


[Vr l+r  1 1 <“y  )|VVyV  I'2 V rv  1 1  V[i Vw  1^ l“yr  1 1  "« 


(25) 
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Figure  11.  Single-Pair  Array 

The  sensitivity  vectors  for  the  accelerometer  can  be  expressed  as 


S.  =(a.  U  +  b.  U  +  c.  U  ) . 
i  i  x  1  y  i  z 


(26) 


The  outputs  are  the  dot  product  of  the  sensitivity  vectors  with  their 
respective  acceleration  vectors. 


0.  =S .  • A. 

ill 


(27) 


Since  the  coordinates  at  position  1  are  parallel  to  the  local  coordinate 
system  and  the  coordinates  at  position  3  are  antiparallel  to  the  local 
coordinate  system, 


S.  =  a  U  +  b  U  +  c  U  ; 
1  lx  1  y  1  z 


(28) 
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(29) 


S  =  -  a.U  -  b.U  -  c„U  ; 
3  3  x  3  y  3  z’ 


and 


22  ••  •  •  , 

0  =  a,(M  +r-r1(o)  +go  )]+  b,[M  +go  oo  r.+2co  r,+G0  r-]+c.[M  +u>  go  r  -go  r  -a>  r  ] 
1  lxlly  z  l  yzxl  zlzlJl  zzxlzlzl 


(30) 


2  2  •  •  • 

0  =a_  [M  -r,  +r,  (go  +oo  ) ]-b_ [M  -go  go  r-  -2go  r  -go  r  ]-c  [M  -go  go  r  r  +oj  r.  ] 

3  3  x  1  1 v  y  z  ' J  3  y  yxl  zl  zlJ  lLz  zxl  yl  z  1 


(31) 


Antiparallelism  has  the  advantage  of  rotational  symmetry  which  simplifies  the 
design  and  fabrication  of  fixtures.  The  outputs  are  converted  to  acceleration 
units  by  applying  the  appropriate  gage  factor  as  follows, 


G. 

l 


(32) 


which  results  in 


r\  r\  ^  ^ 

G.  =  =  d .  [M  -r  -r  ,  ( to  -ho  )]4e,  [M  +w  u  r.+2io  r_-ho  r .  ]+f .  [M  -ho  to  r-w  r  -to  r  ] 

1  P,  l  xll  y  z  l  yyxl  zlzl  l  zzxlzlzl 


y  y 


(33) 


G„  =  ■—  =  d„[M  -r ,  +r ,  ( to  2-ko  2)]-e„[M  -(0  to  r,+2u  r  -<o  -to  <o  r  -hi  r  -to  r  ] 

3  P„  3  x  1  l  y  z  'J  3  y  yxl  zl  zlJ  3  z  zxl  zl  zl 


•  •  • 


(34) 


where 
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and 


di*i; 


ei-kii> 


f  =k 
i  2i> 


or 


drku- 


ei'k2i- 


£i-i; 


etc . 

It  can  be  shown  that  separating  a  specific  vector  component  from  the  output  of 
a  single  accelerometer  is  impossible,  even  if  a  perfectly  uniaxial  device  is 
used.  For  example,  if  d1=1  and  ej  =  fj  »  0, 


Gl'a  1 A1X ~Rx+r  1  ^xUy Ry  ^x  “z Rz  '( Rx+r  1  }  (  “y  ^z  2  )+2\  Rz "2a)z *y+*z Rz \ Ry 


(35) 


Modeling  results  and  previous  measurements  show  that  the  terms  such  as 

^WyWxri^  can  be  sufficiently  large  to  demand  attention.  Also,  the  angular 
components  in  the  other  vector  terms  cannot  be  ignored. 

In  order  to  separate  components,  additional  pairs  of  accelerometers  must 
be  used  and  the  sums  and  differences  of  their  data  are  required  to  separate 
vector  components. 


3A 


(36) 
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v  2  J  y  z  v  2  '  y  x  zlz 


f  ,-f  O 

+  ( — 2 — ')(“z“xr-2<V"“yr* 


G  +G  d  -hi.  ..  ?9  e  -he,. 

— - —  =  ( — z — )(r-r(a>  +  a)  )  +  ( — a)  r+2o)  r+w  r) 
2  v  2  '  y  z  v  2  '  y  x  z  z 


f  i  #  d  -d  ~ 

+  (—2— )(“zV"2V”“yr)  +  (— 2“ )MX 
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(37) 


Upon  examination  of  Eqs .  (36)  and  (37),  the  advantage  of  having  matched 
accelerometers  becomes  obvious.  The  commonality  required  is  in  the  cross-axis 
sensitivity  of  the  two  devices.  It  becomes  further  evident  that  the  lower  the 
maximum  cross-axis  sensitivity,  the  better  the  separation  among  components. 


When 


drd3 


GrG3 


d,  M  +  e.  M  +  f  M 
lx  1  y  1  z 


61  e3 


->•0 


frf3 
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y  x  z  z  9 
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z  x  y  y 


(39) 


when  d,  e  or  f =1 •  Realistically,  what  is  achievable  is  a  pair  of 
accelerometers  with  matched  cross-axis  sensitivity  characteristics  and  maximum 
cross  sensitivities  less  than  1%.  That  is,  =  k^,  1^1  =  k23  anc*  we 

continue  our  analysis  based  on 


kU  <  -01-  ku 


*  k13-  k21 


k23>  and  k21 


<  k 


li* 


Since  it  is  obvious  that  a  single-pair  array  does  not  provide  adequate 
definition  of  single  vector  components,  multiple-pair  arrays  must  be  examined. 

B.  Multiple-Pair  Arrays 

Before  examining  and  analyzing  multiple-pair  arrays,  the  acceleration 
vector  at  a  point  on  the  gun  tube  shown  in  Eq .  (6)  should  be  examined.  To 
separate  M^,  and  ,  one  can  take  advantage  of  the  fact  that  these 
components  do  not  vary  with  the  local  position  vector  component  r.  This  opens 
several  options  to  the  experimentor  in  placing  accelerometer  array  pairs.  The 
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fact  that  the  local  angular  acceleration  components  are  directly  dependent 

on  r  must  be  acknowledged  when  trying  to  separate  rotational  components.  It 

should  be  noted  that  the  components  of  r  are  factorable  in  the  subsequent 

Eqs,  (36)  through  (39),  Eqs.  (36)  and  (37)  show  that  a  minimum  of  three  pairs 

is  required  in  an  array  to  unequivocally  separate  M  ,  M  and  M  ,  where  each 

x  y  z 

pair  consists  of  matched  accelerometers.  The  amount  of  mismatch  of  the  cross¬ 
axis  sensitivity  in  a  given  pair  determines  the  error  included  in  the  data 
analysis.  co^,  and  cannot  be  unequivocally  separated  in  an  array. 

The  multiple-pair  array  can  be  classified  in  three  basic  categories: 
colinear,  coplanar  and  three-dimensionally  distributed. 

1 •  Colinear  Three-Pair  Array .  The  colinear  three-pair  array  is  shown  in 
Figure  12.  Basically, the  origins  of  the  coordinate  systems  of  all  six 
accelerometers  lie  on  the  same  line,  and  their  coordinate  axes  are  parallel  to 

the  local  coordinate  axes  with  those  in  the  upper  position  shown  in  Figure  12 

coparallel  and  those  in  the  lower  position  shown  in  Figure  12  antiparallel. 

In  simpler  terms,  this  array  can  be  visualized  as  a  pair  of  triaxial 
accelerometers  mounted  diametrically  opposite  one  another  on  the  gun  tube  with 
their  respective  axes  opposing  one  another.  For  the  indexing  of  positions, 
the  first  pair  is  at  positions  1  and  3;  the  second  pair  is  at  positions  5  and 

7;  the  third  pair  is  at  positions  9  and  11.  The  coefficients  for  the  array 

are  given  in  Table  5. 


TABLE  5.  COLINEAR  THREE-PAIR  ARRAY  COEFFICIENTS 


Position 

di 

e. 

i 

fi 

ri 

Pi 

i 

1 

kn 

k21 

rl 

P1 

3 

1 

k13 

k23 

"rl 

P3 

5 

k25 

1 

k15 

r5 

P5 

7 

k27 

1 

k17 

-r5 

P7 

9 

k19 

k29 

1 

r9 

P9 

11 

klll 

k2 1 1 

1 

~r9 

P11 
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Figure  12.  Colinear  Three-Pair  Array 

•»  2  2  f  # 

The  vector  components  M  ,  M  ,  M  ,  [r  -r,((D  +o)  ;],[<d  o>  ^+20)  r  ,+U)  r, 

x*  y*  z’  1  1 v  y  z  J,L  y  x  1  z  1  z  1 

[<D  m  r  -2o3  r  -o)  r ,  1  can  be  found  by  taking  the  sums  and  differences  of 
zxlylyl  J  G 

pairs . 
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(40) 


By  substitution,  Eqs.  (40)  become  Eqs.  (41)  when 


(k. ,+k  1 (k  +k0  ) 

J_ijik JK  Im  in 1 


<  .0001 


and  is  considered  to  be  zero  for  measurement  purposes,  and  r  is  considered  to 
be  of  constant  magnitude 
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(k  +k  )+r  +(k  -k  )r  (k  -k  )r  -(k  +k  )r 

Mx-  V2  cr  V2  G3-(  11  13  5  4“  13  1)G5-(  11  V  11  13  5)G7 
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)G7+1/2G9i  *  * 


(k29+k2n 


)r9+(k29~k211)r5'„  r(k29+k211)r9  (k29_k211)r5 

~5  JG5'1  4?^ - 


(41) 


*Not  during  projectile  passage  when  r  and  r  exist 


The  results  shown  in  Eqs.  (41)  are  then  used  to  establish  the  data  processing 
algorithm. 


2*  Coplanar  Three-Pair  Array.  It  is  not  always  physically  possible  to 
arrange  accelerometers  in  a  colinear  array  due  to  other  physical 
constraints.  Therefore,  coplanar  arrays  must  be  considered.  If  an 
accelerometer  is  mounted  on  the  surface  of  the  tube  at  the  end  of  the  position 
vector  r  with  one  principal  axis  coparallel  to  r  and  the  other  principal  axis 
coparallel  to  the  tube  axis,  as  shown  in  Figure  13,  then  the  acceleration  and 
sensitivity  equations  are: 


_  2  2 
V  {[Mx+ri+ri(a)y  +  [Myri“ywz+2Vi+^zri]cos<t»i}ui 

+  {[Hy+r1MyMz+2wzr1-t«zr1l  sin^-lM+^+r^w }u2 


+  {Mz+[«z«yr1-2«xJi-ixr1]co8*i+[«z«xri-2«  r.-w r^sin*.  }u  ; 


(42) 


and 


P. 

l 


dl  U1  +  ei  U2  +  fi  U3 


(43) 


X 


►  Y 


— U  t  —  Ua 

Figure  13.  General  Single-Pair  Array 


The  array  pair  is  completed  by  mounting  an  antiparallel  mate  in  the 
diametrically  opposite  position  i+2 .  The  following  general  equations  for  the 
difference  and  sum  are  generated. 


Gi~Gi+2  si  ;  Ai  Sl+2  1  Ai+2 
2  "  pi  Pi+2 
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S.  •  A.  S.  •  A. 
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Figure  14  shows  an  example  where  the  three-pair  array  is  aligned  with  the 
coefficients  given  in  Table  6. 
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Figure  14.  Coplanar  Three-Pair  Array 


TABLE  6.  COPLANAR  THREE-PAIR  ARRAY 
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By  substitution  and  recasting  the  difference  of  paired  outputs 
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By  closely  examining  the  Eqs .  (46),  it  is  evident  that  M  ,  M  and  M  cannot 
be  separated  without  residual  error  unless  the  corresponding  accelerometer 
pairs  are  matched.  The  use  of  matched  accelerometer  pairs  becomes  a  necessity 
for  a  coplanar  array.  Separation  of  angular  terms  cannot  be  readily  done  with 
real  accelerometers.  These  terms  depend  on  the  maximum  value  of  cross— axis 
sensitivities  being  much  less  than  1%,  say  less  than  .1%. 
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3.  Three-Dimensional  Arrays .  The  difficulties  of  separating  vector 
components  in  three-dimensional  arrays  are  compounded  because  of  the  physical 
restriction  of  being  displaced  from  the  central  axis  of  the  gun  tube.  In 

order  to  circumvent  this  situation,  an  approximate  three-dimensional  technique 

•  • 

is  used  to  determine  and  between  two  three-pair  arrays. 

Consider  a  segment  of  the  gun  tube  as  depicted  in  Figure  15  where  the 

ends  of  the  segment  correspond  to  the  locations  of  two  array  planes  and 

M  ,  M  and  M  are  defined  in  each  plane.  Also,  the  length  L  of  the  segment  is 
X  y  z 

sufficiently  short  to  cause  the  angle  ip,  the  misalignment  angle  between  the 
two  planes,  to  be  insignificant.  Under  these  conditions,  the  rotations  can  be 
formulated  as  follows: 


Figure  15.  Three-Dimensional  Array 
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(47) 


M 


9 


(48) 


Before  Eqs .  (47-48)  can  be  used,  the  length  L  must  be  examined  to  insure 
that  the  previously  stated  assumption  of  sufficient  shortness  is  not 
violated.  The  angle  can  be  expressed  in  terms  of  L  based  on  the  Bernoulli 
beam^  equation 
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El 
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2  2 
E(wZ+l)CT 


,  2  N 
(vz) 


.  dx  rz  n  - 

f'=77=/m  EIdz 
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dz 


8Rp 


E(w2+1)C2 


(49) 


where 


I  =  cross-sectional  moment  of  the  tube  at  z; 

Mq  =  moment  on  the  tube  at  z; 

E  =  modulus  of  elasticity; 
p  =  density; 

M  =  average  transverse  acceleration; 

C  =  caliber  of  the  tube; 
w  =  wall  ratio  of  the  tube; 

zm  =  distance  from  the  primary  array  to  the  muzzle; 

L  =  distance  from  the  primary  array  to  the  secondary  array; 


7 

Ferdinand  L.  Singer,  Strength  of  Materials,  Harper  &  Brothers,  1951. 
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The  constraint  for  the  assumption  is 


8Rp _ 

E(w2+1)C2 


(z 


.001 


(50) 


Eq .  (50)  reduces  to 


E(w  +1)C_  (>001) 
8Rp 


(51) 


Experience  tells  us  that  over  the  range  of  caliber  and  wall  ratios 

encountered  z  should  not  exceed  2.7  calibers  and  the  minimum  L  should  not  be 
m 

less  than  0.5  caliber.  The  condition,  inequality  (51),  has  to  be  verified  for 

each  case.  Eqs.  (48)  then  are  used  in  conjunction  with  the  Eqs.  (41)  and 

•  •  • 

(47),  depending  on  the  types  of  arrays  used,  to  determine  0)  ,  a)  and  u>  • 

x  y  z 

V.  TRANSFORMATION  TO  EARTH  COORDINATES 

In  the  previous  discussion,  the  vector  components  are  described  in  terms 
of  the  local  coordinates  of  the  gun  tube.  That  is,  in  terras  sensed  by  an 
observer  on  the  gun  tube.  However,  it  is  most  convenient  to  describe  the 
motion  of  elements  of  a  gun  system  relative  to  an  arbitrary  Earth  coordinate 
system.  This  type  of  description  allows  the  determination  of  the  relative 
motion  between  various  components  of  the  gun  system. 

A.  Magnitude  of  Rotation  Considerations 

Before  performing  coordinate  transformations,  the  rotational  aspects  of  the 
motion  must  be  examined. 
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0)  dt  +  /  /  0)  dt  dt ; 

x  J  J  x 

o  o  o 
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t  t  t 

ez  =  ez  +  /  coz  dt  +  /  /  co  dt  dt 

o  o  o  o  o 


where 


0^  =  initial  angular  aspect; 

o 


foh  =  initial  angular  velocity; 
o 


(i). 

l 


=  instantaneous  angular  acceleration. 


For  most  experiments  the  gun  system  is  at  rest  at  the  time  t=0.  Consider  Eqs . 

(52)  where  w.  =  0; 

o 


t  t 

9=0  +  /  /  0)  dt  dt ; 

x  x  J  J  x  * 

o  o  o 


t  t 

ey  -  0V  +  /  /  i  dt  dt; 

J  ■'a  r\  r\  * 


t  t 

9  =  0  +  /  /  u  dt  dt. 

Z  Z  J  J  Z 

o  o  o 


(53) 


If  the  absolute  magnitudes  of  0^  0y  and  @z  remain  less  than  0.06  radian, 
then  an  orthogonal  transformation  can  be  used.  0^  ,  0  and  0  must  be 
determined  prior  to  the  measurement  by  other  means.0  °  ° 

If  the  previously  stated  condition  is  violated,  the  transformation  must 
be  performed  through  the  use  of  Euler  angles,  and  the  angular  components  must 

g 

be  redefined. 

5  * 

Henry  L.  Langhaar 3  Energy  Methods  In  Applied  Mechanics ,  John  Wiley  and  Sons3 
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Figure  16.  Orthogonal  Transformation  Angles 
B.  Orthogonal  Transformation 

Figure  16  shows  the  relationship  between  the  local  and  Earth  coordinates  in 

terms  of  the  orthogonal  angles  9,8  and  8  •  Based  on  the  condition  that  the 

x  y  z 

absolute  magnitudes  of  9  ,  8  and  8  are  less  than  0.06  radian,  the 

x  y  z 

transformation  matrix  is  shown  in  Table  7  where  f(U  ,U  ,U  )  is  the  local 
a  a  a  x  y  z 

vector,  and  F(U  ,U  ,U  )  is  the  Earth  vector. 

*  x  y  z 

TABLE  7.  ORTHOGONAL  TRANSFORMATION  COEFFICIENTS 
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(54) 


M  =  M 
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__  ^ 

U  +M  U  +M  U 
x  y  z  z 


(M  -0  M  +0  M  )U 
x  z  y  y  z'  x 


(M  +0  M  -0  M  )U  +  (M  -0  M  +0  M  )U 
yzxxzy  zyxxy  z 


The  components  of  motion  in  the  Earth  coordinate  system  become 
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(M  -  0  M  4-  0  M  )  ; 
x  z  y  y  z  * 


M  -(M+0M-0M); 
y  y  z  x  x  z'* 


M  '  =  (M  -  0  M  +  0  M  )  . 
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(55) 


0  =0  =0  =  0:  when  t  =  0. 

x  y  z  5 

o  o  o 

0  =  A;  0  =  B;  0  =  C;  when  t  =  0. 

x  y  z 

o  o  o 

C.  Transformation  By  Euler  Angles 

Figure  17  shows  two  arbitrary  coordinate  systems  rotated  with  respect  to 
each  other  by  the  angles  0,  <j>  and  Table  8  shows  the  direction  cosines 
between  the  unit  vector  of  the  coordinate  systems.  This  matrix  is  valid  as 
long  as  0*0.  Numerically,  the  transformation  breaks  down  when  0  is  in  the 
neighborhood  of  zero.  By  judicious  selection  of  the  coordinates,  0  can  be 
made  to  be  always  in  the  neighborhood  of  tt/2. 
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TABLE  8.  GENERAL  EULER  ANGLE  TRANSFORMATION  COEFFICIENTS 


n 

c 

X 

sin(c())sin(iJO 

-cos  (9)cos  (c())cos(\|;) 

sin(c())cos(^) 

+cos  (9)cos  (c())sin(^) 

sin(9)cos(c()) 

y 

-cos(c())sin(\|;) 

-cos  (  0)sin(c{))cos  (\|>) 

-cos(c())cos(^) 

+cos(9)sin(c())sin(iJ;) 

sin(9)sin(c()) 

z 

sin(9)cos(\J;) 

-sin(9)sin(\{0 

cos(9) 

Y  Z 


Figure  17.  General  Euler  Angles 


Figure  18  shows  the  coordinate  systems  for  the  desired  transformation.  A  new 
table  can  be  constructed  by  making  the  following  substitutions. 
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(56) 


Figure  18.  Euler  Angle  Transformation  of  local  and  Earth  Coordinates 

Table  9  shows  the  new  transform  matrix  after  the  Eq.  (55)  has  been 
substituted  in  the  matrix  in  Table  8. 
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D.  Six-Degree-of -Freedom  Measurement 


On  a  gun  tube,  the  six-degree-of-f reedom  measurement  can  be  accomplished  by 
using  two  colinear  three-pair  arrays. 

1.  Limitations .  This  measurement  is  limited  to  positions  close  to  the 
muzzle . 

2.  Requirements .  This  measurement  requires  that  the  three-pair  arrays 
be  colinear.  The  positions  of  the  two  arrays  conform  to  Eq.  (51). 

E.  Linear  Three-Degree-of -Freedom  Measurements 


The  linear  vectorial  components  of  tube  motion  can  be  determined  by  using 
either  colinear  or  coplanar  arrays. 

1.  Limitations .  Three-degree-of-f reedom  measurements  do  not  allow 
transformation  to  Earth  coordinates.  Coplanar  arrays  can  be  used  as  well  as 
colinear  arrays, but  must  have  matched  pairs  of  accelerometers. 

2.  Requirements .  Colinear  arrays  must  have  matched  pairs  of 
accelerometers.  Positions  and  orientations  must  be  well  controlled  and 
recorded . 

F •  Data  Processing 


As  can  be  seen  from  the  equations  and  conditions  which  must  be  used  to  make 
acceleration  measurements,  the  subsequent  processing  should  not  be  attempted 
without  extensive  automated  facilities  to  digitize  and  numerically  process 
massive  amounts  of  data. 

One  of  the  most  significant  aspects  of  the  analysis  results  is  that  data 
processing  system  bandwidth  must  accommodate  the  second  harmonic  of  the 
highest  mode  frequency  of  the  structural  response.  This  means  that 
digitization  rates  must  be  twenty  times  the  highest  response  frequency  of 
interest  in  order  to  adequately  depict  the  six-degree-of-f reedom  motion  of  gun 
and  projectile  systems.  This  requirement  also  applies  to  the  data  rate  and 
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time  step  increment  of  numerical  simulations  of  gun  and  projectile  system 
response . 

G.  Comparison  of  Measurement  and  Theoretical  Predictions 

The  primary  result  of  the  preceding  analysis  is  that  simulation  outputs  must 
be  given  in  terms  of  the  components  that  can  exist  in  the  measurand,  not  in 
the  components  historically  used  for  mathematical  convenience j  that  is 


R-KjOojR-KjOojR-R 
x  xyy  x  z  z  x 


(u)  -kd  )+2o)  R  -2o)  R  *Ho  R  -a)  R  =M 
x  y  yz  zyyzzyx 


a  viable  measurand  where 


R 

x 


is  not. 
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